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Abstract 

We compare the behavior of different lattice Dirac operators in 
(~| , gauge backgrounds which are lattice discretizations of a classical in- 

stanton. In particular we analyze the standard Wilson operator, a chi- 
rally improved Dirac operator and the overlap operators constructed 
■ from these two operators. We discuss the flow of real eigenvalues as a 

function of the instanton size. An analysis of the eigenvectors shows 
that overlap fermions with the Wilson operator as input operator have 
difficulties with reproducing the continuum zero mode already for mod- 
erately small instantons. This problem is greatly reduced when using 
the chirally improved operator for the overlap projection. 
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Introductory remarks: Recently it was realized that the Ginsparg- Wilson 
equation [P is crucial for implementing chiral symmetry on the lattice. Cur- 
rently three types of exact solutions are known: The overlap operator 
perfect actions Q and domain wall fermions Furthermore a systematic 
expansion of a solution of the Ginsparg- Wilson equation was developed and 
tested in [|, |]. These new lattice Dirac operators based on the Ginsparg- 
Wilson equation were recently used to analyze relevant excitations of the 
QCD vacuum which affect the Dirac operator ^, ^, 10|, in particular the 



local chirality variable proposed in |11] was studied in detail. 

These studies of the excitations affecting the lattice Dirac operator are 
motivated by the instanton picture of chiral symmetry breaking (see [^] 
for extended reviews). A single instanton or anti-instanton produces a zero 
eigenvalue of the Dirac operator. An interacting pair of an instanton and 
an anti-instanton leads to a complex conjugate pair of small eigenvalues 
instead of two zero eigenvalues. In instanton models the QCD vacuum is 
pictured as a fluid of interacting instantons and anti-instantons which lead 
to an accumulation of small eigenvalues near the origin. Since the density 
of eigenvalues at the origin is related to the chiral condensate through the 
Banks-Casher relation the fluid of instantons and anti-instantons leads 
to a breaking of chiral symmetry. The lattice studies 0-|lC],|14] tried to 
prove or refute this picture of interacting instantons and anti-instantons. To 
be more specific, various observables built from the eigenvectors of the Dirac 
operator were studied for background gauge fields generated by simulations 
in the quenched approximation. 

A good method for testing properties of different lattice Dirac operators 
is to study them in smooth instanton backgrounds. In this letter we report 
on such a study comparing the overlap Dirac operator based on the Wilson 
operator Q, the standard Wilson operator, a recently proposed approxi- 
mate solution of the Ginsparg- Wilson equation |^, § which we will refer to 
as the chirally improved Dirac operator and the overlap operator with the 
chirally improved operator as input operator. We analyze different prop- 
erties of the eigenvectors and eigenvalues for these Dirac operators using a 
lattice discretization of instantons. Another study including smooth instan- 
ton backgrounds was reported in , which investigated the number of zero 
modes and the topological charge as a function of instanton radius. 

The goal of our analysis is twofold: Firstly it serves to better understand 
the results of the above mentioned studies of relevant excitations in the QCD 
vacuum as seen by the lattice Dirac operator. When the instanton is large 
compared with the lattice spacing, all operators give good results. But when 
we consider smaller instanton radii we find that the overlap operator with 
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Wilson input operator has difficulties with reproducing the continuum zero 
mode in an instanton background even for moderately small instantons. We 
show that this problem is greatly reduced when using the chirally improved 
operator for the overlap projection. 

The second goal of this study is to achieve a better understanding of 
technical aspects of the overlap projection. The sensitivity to defects should 
be understood, in particular if one uses an approximate solution of the 
Ginsparg- Wilson equation as a starting point for the overlap projection as 



proposed in |16| 



Technicalities: We analyze the overlap Dirac operator, the standard Wil- 
son Dirac operator and the chirally improved Dirac operator. The latter 
has been described in detail in Here, since the discretized instanton 
configurations are smooth, we use the coefficients for the free case as listed 
in the appendix of [^. 

The overlap operator is given by 

P.. = 1- , (1) 
[l-Dl){l-Do) 

where Dq is any reasonable lattice version of the Dirac operator. Here we 
use for Dq both the standard Wilson Dirac operator leading to the Wilson 
overlap operator as well as the chirally improved operator leading to the 
chirally improved overlap operator. We construct the inverse square root 
using Chebychev approximation following the approach discussed in [17, 



A method for putting an instanton on the lattice has been proposed in 
Ref. |19|. Here we use a slightly modified version of this procedure. In the 



regular gauge, the continuum instanton potential is written as 



^m(^) ~ 2(x^ + p^) ^^^^^^ SyS^)xu, (2) 



and in the singular gauge we have 

^m(^) = 2a;2(3;2 + p^) ^^^'^'^ ~ SyS^) Xu , (3) 

where S4 = S4 = 1, = — = \aj {j = 1, 2, 3) and p is the "radius" of the 
instanton. The corresponding expressions for an anti-instanton are obtained 
by exchanging sx and sx. 

The first step of our procedure is a coordinate transformation x^ — > 
2/^ which maps the real line onto the interval (0,L), where L is identified 
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with the length of our lattice (measured in units of the lattice spacing and 
assumed to be even). For this transformation we take = /(y^) with 

f{y) = L^[{L-y)-^-y-^] . (4) 

Acting with this transformation on the (anti-)instanton potential either in 
the regular or in the singular gauge we obtain the corresponding potentials 
on the four-torus (0, L)^. The center of the instanton, xi = • ■ ■ = X4 = 0, is 
mapped onto yi = • • • = 7/4 = L/2. Note that the instanton is "squeezed" 
by this procedure so that the radius R of the potential on the four-torus is 
related to the radius p of its infinite- volume precursor hy p = f{R + L/2). 

In the second step we divide our lattice into an "inner part" around 
yi = ■ ■ ■ = 1/4 = L/2 and a complementary "outer part". According to the 
procedure of Ref. |jl9| we work with the potential in the regular (singular) 
gauge in the inner (outer) part using the gauge transformation connecting 
the two gauges to glue both potentials together. The third step consists in 
computing the gauge links from the potential. This is easily done analyti- 
cally. 

The resulting SU(2) link variables are finally embedded in SU(3) in the 
most trivial way, namely as 2 x 2 blocks in the upper left-hand corner of the 
SU(3) matrices. 

We work on lattices with size 16^ and 12^. For the fermions we use 
periodic boundary conditions in space direction and anti-periodic boundary 
conditions in time direction. The numerical computations of eigenvalues and 



eigenvectors are done with the implicitly restarted Arnoldi method |2C]. 



The flow of real eigenvalues: It is known that in an instanton back- 
ground an exact solution of the Ginsparg- Wilson equation has exact zero 
eigenvalues j^, as does the continuum Dirac operator. Currently only the 
overlap operator shows this property. A practical implementation of per- 
fect actions requires a finite parametrization of the lattice Dirac operator 
and the coefficients of the parametrization are determined from renormaliza- 
tion group transformations. Since an exact solution of the Ginsparg- Wilson 



equation is necessarily non-ultralocal |21|, any practical implementation of 



the perfect action will be an ultralocal approximation of a solution of the 
Ginsparg- Wilson equation. Thus the fixed point operator as well as the chi- 
rally improved operator will only have approximate zero modes. Both these 
Dirac operators obey 75-hermiticity, i.e. 75-D75 = . This implies |22] that 



eigenvectors ip oi D with eigenvalues A have V'WsV' = unless A is real. This 
has to be compared with the property that V'^TsV' = unless A is zero, which 



3 




Figure 1: The dependence of the position x of the real eigenvalue (zero 
mode) on the radius R (in lattice units) of the underlying instanton. We 
show our results for the Wilson overlap operator (diamonds), the chirally 
improved operator (circles) and the Wilson operator (squares). The data 
were computed on 16^ lattices. 



holds for eigenvectors of an exact solution of the Ginsparg- Wilson equation 
or for the eigenmodes of the continuum Dirac operator. It thus follows that 
for perfect actions, for the chirally improved Dirac operator and also for Wil- 
son fermions only eigenvectors with real eigenvalues are possible candidates 
for topological modes. 

An interesting question is, how well different operators manage to project 
the real mode into the origin when the underlying gauge field changes. In 
Fig. H] we show the position x of the real eigenvalue as a function of the 
radius R (in lattice units) of the underlying instanton configuration. This 
study was done on lattices of size 16^. We display data for the Wilson 
overlap operator, the chirally improved operator and the standard Wilson 
operator for identical gauge configurations. The behavior of the zero modes 
of the Wilson overlap operator is of course trivial and it serves only as a 
reference line. For the Wilson operator we find a very strong dependence of 
the real mode on the radius of the instanton. Already for large instantons 



4 



the corresponding eigenvalue is shifted to relatively large real values and this 
shift increases as the radius of the instanton shrinks further. The chirally 
improved Dirac operator is considerably less sensitive to the radius of the 
instanton. It starts to deviate from only for radii below 2.5 lattice units. 

Our analysis sheds light on a potential problem of the overlap projec- 
tion: Whenever the background configurations contain defects, i.e. gauge 
configurations with very small excitations carrying topological charge, the 
overlap projection becomes numerically expensive and low eigenvalues of 
(1 — -Do)(l — -Do) have to be projected out before the square root in (||) can 
be evaluated numerically. The underlying mechanism is nicely illustrated in 
our Fig. |l|: For small instanton radius the real eigenvalue of the Wilson op- 
erator Dq used in the overlap projection (|l]) comes close to the center of the 
projection (1 in the complex plane) causing the inverse square root to blow 
up which spoils the numerical evaluation of the overlap operator. Comparing 
the Wilson curve with the curve for the chirally improved operator shows 
that when using already an approximate solution of the Ginsparg- Wilson 
equation as Dq in the overlap projection the problem with defects is 
milder. 




Figure 2: Schematic picture for the movement of the physical real mode and 
its doubler partner as a single instanton is destroyed. The circle represents 
the Ginsparg- Wilson circle in the complex A plane, and the two other curves 
are the schematic trajectories of the two eigenvalues. 



Fig. H] also illustrates how topological modes are treated by solutions 
and approximate solutions of the Ginsparg- Wilson equation. For an exact 
solution of the Ginsparg- Wilson equation it is known that the spectrum 
depends discontinuously on the underlying gauge field. Since (1) the total 
number of eigenvalues is even, (2) all eigenvalues have to lie on the Ginsparg- 
Wilson circle (the circle with radius 1 and center 1 in the complex plane) 
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and (3) all modes which are not real come in complex conjugate pairs, a sin- 
gle eigenvalue has to vanish discontinuously as the underlying gauge field 
is deformed from topological sector 1 to sector 0. Since an ultralocal ap- 
proximation of a solution of the Ginsparg- Wilson equation cannot show such 
discontinuous behavior the change of the topological sector of the underlying 
gauge field has to manifest itself differently: As the instanton shrinks, the 
real mode starts to travel into the interior of the Ginsparg- Wilson circle (for 
an approximate solution of the Ginsparg- Wilson equation the eigenvalues 
are not confined to the circle) where it meets a partner from the doubler 
branch of the spectrum. When they meet on the real axis they can contin- 
uously form a complex conjugate pair and travel back to the outside of the 
circle. A schematic picture of this behavior is given in Fig. |2[ 

The better a Dirac operator approximates a Ginsparg- Wilson fermion, 
the faster the eigenvalue moves through the center of the circle. We remark 
that we have also seen the behavior of Fig. |2|, in which a partner from the 
doubler sector meets the displaced real mode to form a complex conjugate 
pair, in a numerical study of random matrices with the same symmetries as 
the lattice Dirac operator. 

Localization properties of the eigenvectors: In order to further an- 
alyze the behavior of different lattice Dirac operators in instanton back- 
grounds we now study properties of their eigenvectors. In the background 
of an instanton field the continuum Dirac operator has a zero mode tpo (see 
e.g. fl^). It is localized at the same position as the underlying instanton. 
A gauge invariant density p{x) which inherits this localization is obtained 
by summing over the color and Dirac indices c and a, 

2E? 

a,c -K yn ^ X ) 

Due to the normalization of the zero modes we have / d'^xp{x) = 1. A 
measure of the localization of ip is given by the inverse participation ratio, 

I = J (fxp{xf = (57r2i?4)-i. (6) 

As an alternative measure of locality ||l^ studies the self-correlation of the 
local chiral density. For different radii of our lattice instantons we computed 
the inverse participation ratio of the corresponding zero mode. In the plot 
on the left-hand side of Fig. ^ we show our results for lattice size 16^. We 
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Figure 3: The inverse participation ratio of the zero mode as a function of 
the radius R of the underlying instanton. On the left-hand side we show 
our results on IG'^ lattices for the Wilson overlap operator (diamonds), the 
chirally improved operator (circles) and the Wilson operator (squares). The 
right-hand side plot shows a comparison of the Wilson overlap operator 
(diamonds) and the chirally improved overlap operator (asterisks) on 12^ 
lattices. The dashed lines represent the continuum formula. We rescale / 
by the volume, i.e. we plot / x 16^, respectively / x 12'*. 



normalized the inverse participation ratio / by the volume^, i.e. we plot / x 
16'*. The symbols give the numerical results while the dashed line represents 
the continuum formula from Eq. (^). From the plot is obvious that for 
our largest value of the instanton radius, R = 6, the lattice results are 
slightly above the continuum value due to finite size effects. For smaller R, 
the numerical results for the two ultralocal lattice Dirac operators, i.e. the 
Wilson Dirac operator (squares) and the chirally improved Dirac operator 
(circles) follow the continuum curve down to relatively small values of R. 
For the chirally improved operator the agreement with the continuum result 
holds down to ii = 1 and only for R = 0.75, as the instanton begins to 
"fall through the lattice" we find a considerable deviation. For the Wilson 
operator, due to the move of the real mode into the interior of the eigenvalue 
distribution we could not obtain data for R < 2. However, down to i? = 2 

^With this normalization the inverse participation ratio is a quantity widely used in 
solid state physics. 
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Figure 4: The maximum Pmax of p{x) as a function of the radius R of the 
underlying instanton. In the plot on the left-hand side we show our results 
on 16^ lattices for the Wilson overlap operator (diamonds), the chirally 
improved operator (circles) and the Wilson operator (squares). The right- 
hand side plot provides a comparison of the Wilson overlap operator and 
the chirally improved overlap operator on 12^ lattices. 



the Wilson operator also follows the continuum curve quite well. 

The situation is different for the Wilson overlap operator. Already at 
R = 3.5 the zero mode of the Wilson overlap operator has a value of the 
inverse participation ratio which is visibly different from the continuum re- 
sult, and the error quickly increases as R is decreased (note the logarithmic 
scale). At -R = 2.5 the Wilson overlap result amounts to only 60% of the 
continuum formula and at i? = 1.5 only about 30% remain. 

In the plot on the right-hand side of Fig. |^ we compare the inverse 
participation ratio of the Wilson overlap operator (diamonds) to the inverse 
participation ratio of the zero mode of the chirally improved overlap operator 
(asterisks). It is obvious that with the chirally improved overlap operator 
the results for small instantons are considerably closer to the continuum 
line. 

A second characteristic quantity of the zero mode is the maximum Pmax 
of p{x), which in the continuum is given by Pmax = p(0) = 2/{tt'^R'^). In 
Fig. ^ we plot our data for this quantity as a function of the instanton radius 
R. 
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Again the left-hand side plot shows the results for the chirally improved 
operator, the Wilson operator and the Wilson overlap operator, while the 
right-hand side compares the Wilson overlap operator with the chirally im- 
proved overlap operator. The overall picture is similar to the results for 
the inverse participation ratio. The chirally improved operator gives the 
best results, while the Wilson overlap operator has again problems with re- 
producing the continuum results for smaller instantons. One finds that the 
amount of the deviation of the Wilson overlap result from the continuum 
formula is already 50% at i? = 2.5 and increases further for smaller R. The 
situation is improved when using the chirally improved overlap operator. 

We remark that we performed the same analysis with two changes of 
our setting: (1) Instead of using 1 as the center for the overlap projection 
we also used 1 + s, with s = 0.1, s = 0.2 and s = 0.5. Such an adjustment 



of the center of the projection is known |17] to optimize the localization 
properties of the Wilson overlap operator. We found that a variation of 
s amounts to only small changes of / and Pmax for the zero modes of the 
Wilson overlap operator. (2) We also used a different discretization of the 
continuum instanton. We placed the instanton at the center of a hypercube. 
This allowed us to shrink the "inner part" of the gauge potential such that 
it only consisted of the interior of this hypercube, and we could then use the 
potential in the singular gauge on the whole lattice. Also this modification 
does not change the picture we obtained and shows that the results are not 
very sensitive to the details of the discretization of the instanton. 

To summarize, we find that the Wilson overlap operator has significant 
problems with reproducing the continuum zero mode for instantons with 
R < 2.5, while the other operators, i.e. the chirally improved operator, the 
Wilson Dirac operator and the chirally improved overlap operator do not 
show such a large deviation from the continuum result. 

Why does the chirally improved overlap operator reproduce the contin- 
uum results better than the Wilson overlap Dirac operator? Possibly this 
is connected with the fact that the input Dirac operator in the chirally im- 
proved case is already a much better approximation to a Ginsparg- Wilson 
operator than the Wilson operator. 

For a typical simulation with a ~ 0.1 fm, our results for the Wilson over- 
lap operator imply that structures smaller than ~ 0.3 fm will probably not 
be resolved properly. We expect that the chirally improved operator with or 
without additional overlap projection fares better for such small structures. 
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